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ABSTRACT 


The  theoretical  basis  for  the  CAFA  program  is  discussed. 
An  approximate  technique  evolving  from  the  theory  is  applied 
to  the  analysis  of  the  current-voltage  characteristic  of  a 
hypothetical  diode,  with  good  results.  A  printout  of  the  re¬ 
sultant  program  and  data  is  included. 
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INTRODUCTION 


The  CAFA  {Computer-Aided  Fit  and  Analysis;  pronounced 
”ca£6")  program  was  developed  for  fitting  sihooth  curves  to 
experimental  current-voltage  and  similar  data  and  for  using 
the  curves  obtained  to  analyze  the  data.  The  original  goal 
was  to  allow  the  slope  of  such  curves  to  be  obtained  at  any 
point  in  order  to  assist  in  determining  the  current  mechanisms 
existing  in  solid-state  devices .  The  program  was  found  to  be 
useful  for  several  other  purposes.  One,  in  particular,  allowed 
interpolation  between  data  points.  The  preliminary  results 
obtained  have  been  encouraging.  The  program  was  used  exten¬ 
sively  to  analyze  the  data  discussed  in  Reference  1. 

The  subroutine  which  makes  the  program  possible  is  the 
SMOOTH  subroutine  originally  developed  by  Reinsch  (1967) , 
adapted  for  Fortran  by  R.  E.  Jones  of  Sandia  Laboratories,  and 
modified  by  the  author  to  include  interpolation.  This  sub¬ 
routine  fits  a  series  of  spline  (cubic)  functions  to  the  data 
points  and  smooths  the  transitions  between  functions  by  re¬ 
quiring  continuity  of  the  first  and  second  derivatives  within 
a  certain  error  chosen  by  the  user.  The  first  and  second  de¬ 
rivatives  are  available  as  printout  in  addition  to  the  fitted 
data  points  and  the  coefficients  of  each  cubic  equation.  The 
latter  permitted  interpolation  between  data  points.  Examina¬ 
tion  of  the  program.  Fig.  2,  reveals  the  use  of  these  features. 
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Several  general  approaches  to  determining  current  mecha¬ 
nisms  from  the  CAFA  program  are  discussed  below,  along  with 
the  specialized  version  used  in  this  study. 


2 


THEORETICAL  BACKGROUND 


Suppose  the  current  in  a  diode,  I,  is  given  by 

I  »  IQ  exp (6V/m)  ,  (1) 

where  I  is  constant,  B  =  q/kT,  V  is  the  applied  voltage,  and 
m  has  a  constant  value.*  In  this  case,  the  current  in  the 
diode  is  described  by  one  mechanism  and  the  slope  of  the 
ln(I)  vs.  V  curve  is 

91n(I)  _  31 

— 5v - 5T  ’  {2) 

Since  everything  in  Eq.  (2)  is  known  except  m,  m  can  be  deter¬ 
mined  and  the  current  mechanism  can  often  (but  not  always)  be 
identified.  Unfortunately,  total  diode  currents  given  by  Eq. 
(1)  usually  do  not  exist  in  practice.  The  current  in  real 
diodes  is  more  often  given  by  an  expression  of  the  form 

I  =  l  IQi  exp(BV/mi)  +  l  I  exp(V/<J>.)  ,  (3) 

i  j  J  J 

where  describes  the  i'"  temperature-dependent  mechanism 

(e.g.,  diffusion,  space-charge  region  recombination,  surface) 

and  <J>j  describes  the  j*"*1  nontemperature-dependent  mechanism 

(e.g.,  tunneling).  Even  Eq.  (3)  does  not  describe  the  most 

general  current  form  because  it  neglects  nonlinear  effects, 

such  as  interactions  between  current  components,  plus  it 

*To  be  perfectly  general,  "current"  should  be  replaced  by 
"flux"  and  "voltage"  should  be  replaced  by  "force." 
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the  junc- 


assumes  all  the  I  's  are  constant  and  that  V  =  V., 
o  3 

tion  voltage.  Nevertheless,  Eq.  (2)  is  often  an  excellent 
approximation.  For  real  diodes,  one  often  does  not  know  what 
the  current  mechanisms  are  but  could  deduce  the  mechanisms  if 
the  m.'s  in  Eq.  (3)  were  known.  If,  however,  the  diode  cur¬ 
rent  consists  of  only  two  components  and  is  of  the  form 

1  =  X01  eXP(S7>  +  I02  '  (4) 


then 


ain(I)  _  B,l_  1_, 

S\7  I  Vm1  m2; 


(5) 


and  the  two  values  of  m  cannot  be  determined  easily  unless 
each  current  type  has  a  clearly  defined  region  of  dominance 
and  the  experimenter  has  data  covering  those  regions.  Even 
the  simple  case  of  Eq.  (5)  is  not  usually  seen  in  practical 
diodes.  To  further  complicate  matters,  the  value  of  m  describ¬ 
ing  one  mechanism  may  vary;  as  an  example  2  <_  m  <  ®,  depending 

on  injection  level,  for  donor-acceptor  pair  (DAP)  recombination 

2 

in  the  space-charge  region.  A  method  for  finding  the  value 
of  m  at  any  point,  given  a  current-voltage  curve,  is  therefore 
desirable.  A  graphical  approach  suffers  on  several  counts: 

It  is  not  accurate  enough  unless  the  In (I)  vs.  V  curve  is  quite 
linear,  as  will  be  seen  below;  and  it  is  extremely  tedious  and 
time  consuming  to  obtain  m  at  many  points.  The  numerical  ap¬ 
proach  discussed  here  is  quick  and  gives  good  results  for  the 
special  cases  discussed. 
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Given  an  experimental  current-voltage  curve  in  which  the 
current  is  presumed  to  arise  from  the  mechanisms  described  in 
Eq.  (3) ,  we  assume  that  the  current  at  any  point  on  the  curve 
can  be  written  as 


1  =  Io  exp  (mfrr)  '  (6) 

Thus,  our  basic  assumption  is  that  the  experimentally-determined 
current  given  by  Eq.  (6)  is  equivalent  to  the  theoretical  cur¬ 
rent  given  by  Eq.  (3) .  The  parameter  that  allows  Eq.  (6)  to 
describe  the  current  at  any  point  is,  of  course,  m(I) .  If  m(I) 
is  constant,  Eq.  (6)  reduces  to  Eq.  (1) .  The  slope  obtained 
from  Eq.  (6)  is 


31  3  /  V  \ 

w  9  w  (sm;  • 


But  m(I)  =  m(V)  implicitly,  since  I  «  V.  That  is. 


3  V  _  „  3  1  .1  3V 

3v  mTTT  W  mTT 7  m(I)  W  ' 


or 


3  V 

3v  mTTT 


1  3m(I)  3 1 ~|  .  1 

HI?  -n-  5vJ  mnr  • 


(8) 


The  dependence  of  m  on  I  will  be  understood  unless  otherwise 
stated. 

A  differential  equation  can  be  obtained  that  would  solve 
for  m  in  closed  form  if  a  closed  form  solution  exists  for  the 
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differential  equation.  By  combining  Eqs.  (8)  and  (7)  we  obtain 


31 

W  “ 


fll 


m  + 


8 IV 
m 


3m 

7T 


(9) 


By  rearranging  Eq.  (9)  we  obtain,  since  3m/3I  =  dm/di  (m  is  a 
function  only  of  the  current  at  constant  temperature) , 

dm  ±  1  „2  1  _  A  ,1M 

dl  +  FlV  m  V3I/3V  "  0  '  (10) 

This  is  a  very  nonlinear  d.e.  of  the  form 


p.  +  JL.  -  = 

dx  axz  bz 


0  , 


(11) 


where  x  =  x(z)  is  known  and  b  (=  3I/3V)  is  known.  The  d.e. 
might  be  solvable  using  numerical  techniques. 

We  shall  briefly  discuss  the  errors  which  result  when  the 
nonlinearity  in  Eq.  (7)  is  assumed  negligible.  Considering  the 
derivative  on  the  right-hand  side, 


3  ,V.  V  3m  1 

V  7  3V  m  ' 

m 


(12) 


which  is 


3 

Tv 


(V)  «  1 

m  m 


if  3m/3V  is  negligible.  This  requires 


(13) 


6 


V  3m 

T  W  ' 
m 


so  that 


m  >>  -  V  ^  .  (14) 

For  m  x  2,  V  ~  1  volt, 

<<  2  units/volt.  (15) 

The  rate  of  change  of  m  with  V  must  be  very  small  for  this 
condition  not  to  be  violated.  Hence,  the  approximation  that 
3m/3V  is  negligible  is  often  not  valid.  This  approximation 
is  numerically  identical  to  that  obtained  using  an  incremental 
approach  which  approximates  the  curve  between  a  series  of 
closely  adjacent  data  points  by  a  set  of  straight  lines,  and 
is  approximately  equivalent  to  the  results  which  would  be  ob¬ 
tained  using  a  graphical  technique.  Therefore,  unless  a  In (I) 
vs.  V  curve  is  very  linear,  a  graphical  or  incremental  approach, 
or  an  approach  which  neglects  3m/3V  will  not  suffice  for  deter¬ 
mining  the  value  of  m  from  an  experimental  curve. 

As  a  special  case  of  the  relationship  between  Eqs.  (3) 
and  (6)  we  write  the  relation  as 

I  =  l0  exp(^)  =  I01  exp(^)  +  IQ2  exp(££).  (16) 

The  two  terns  on  the  right-hand  side  of  Eq.  (16)  may  be  con¬ 
sidered  as  two  distinct  components  or  as  one  distinct  component 


7 


plus  a  sum  of  other  components,  so  that  the  relation  is  not 
necessarily  restricted  to  only  two  distinct  components.  This 
will  be  called  the  "two-process"  model.  We  shall  assume  Eq. 
(16)  holds  and  determine  the  effective  m  vs.  V  for  various 
ratios  R, 


R  =  101/'I02  * 


(17) 


R  will  be  constant  or  nearly  so  for  many  situations.  Using 
Eq.  (17)  we  can  rewrite  Eq.  (16)  as 


I  =  I 


02 


[R  exp (~)  +  •xp(lj)]  • 


(18) 


The  derivative  is 


31  _  T 
3V  02 


[k  exp®  +  i-  exp®]  . 
Im^  c  m^  m2  m2  J 


(19) 


Define 


f  '  3I/3V  ; 

combining  Eqs.  (18)  and  (19)  gives 

f  =  - 

8 

Equation  (21)  can  be  solved  for  R: 


R  exp(^)  +  exp(j-p) 
m^  mj 


'5-  *  -i  exp  (j-—)' 

m  ^  iti  ^  1x^2  1^2 


(20) 


(21) 


8 


r 


i 

j 

i 


6f 

m 


S±  -  1 


“2 

F"-i 

m. 


exp 


[6vt  -  k] 


(22) 


Since  R  is  constant  with  V,  we  must  have 


3R 

W 


0  . 


(23) 


Performing  the  indicated  operations  on  Eq.  (22) ,  realizing  that 
exp  |  tsv  - |  %  0 ,  and  rearranging  gives 


BV(|-  "  sH 
m  2 


3f 

w 


M  -  i)  {M  -  i)  . 


m 


m. 


(24) 


This  is  the  criterion  for  R  to  be  constant.  One  way  to  use 
Eq.  (24)  would  be  to  find  f (V)  (from  the  results  of  curve¬ 
fitting  to  the  I-V  characteristic  using  the  SMOOTH  subroutine, 
since  this  gives  3I/3V  also),  curve  fit  to  the  points  f,  find 
3f/3V  from  the  curve  fit  and  plot  the  right-hand  side  and 
3f/3V  vs.  Bf  on  the  same  curve,  using  m1  and  m2  as  parameters; 
one  could  inspect  the  results  to  find  integer  values  of  m^  and 
m2  such  that  the  curves  intersect.  Another  way  would  be  to 
find  a  second  expression  describing  3f/3V.  By  differentiating 
Eq.  (20)  directly,  we  obtain 


3f 

37 


l 


32I 


(3I/3V)  3V 


(25) 


Since  the  second  derivatives  can  also  be  obtained  from  the 
curve  fit  of  I  vs.  V,  we  have  enough  information  to  find  3f/3V; 
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that  is,  if  the  second  derivatives  are  reasonably  well  be¬ 
haved.  A  third  method  is  to  rearrange  Eg.  (24) ,  if  3f/3V  can 
be  found,  to  obtain 

m.  +  m~ 

a  £  _  -L  ^ 


Both  sides  of  Eq.  (26)  can  be  plotted  vs.  V  with  m1  and  m2  as 
parameters .  The  solutions  would  be  obtained  at  the  curve  inter¬ 
sections.  For  the  special  case  when  3f/3V  =  0,  Eq.  (26)  re¬ 
duces  to 


(m. 


m2)' 


-  4m1m2(l 


TTZ) 

3V' 


(26) 


6f  =  m^  or  m2  .  (27) 

Thus,  one  of  the  m's  can  be  found  if  3f/3V  *  0  somewhere.  If 
3  f/3V  =  1, 

6f  *  m^  +  m2  or  0  .  (28) 

If  3f/3V  *  1  and  8f  ^  0,  then 

m2  =  6f  -  m1 


if  m1  was  found  at  a  point  where  3f/3V  =  0.  Still  another  ap¬ 
proach  uses  f  in  the  differential  equation  (10) .  Using  the 
definition  of  f,  Eq.  (10)  becomes 


V 


dm 

3v  " 


m 


(29) 


If  $f  is  a  determinable  function,  this  can  be  solved  to  get  a 
family  of  curves  for  various  V's.  Even  if  f  is  not  a 
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determinable  function,  the  expression  (26)  could  be  substituted 
into  the  d.e.  and  then  one  would  have  a  d.e.  in  terms  of  m,  nu  , 
and  m2,  so  that  m  could  be  predicted  given  m^  and  m2  (a  trial- 
and-error  approach) .  For  the  special  case  when  3f/3V  *  1,  the 
d.e.  becomes 


dm 

3v 


(30) 


where  we  have  used  Eq.  (28)  with  (5f  f  0.  This  can  be  solved, 
yielding 


where  m  is  a  known  value  of  m  at  V  =  V 
o  o 

equation  may  also  be  solved  graphically. 


This  transcendental 


A  somewhat  simpler  and  less  tedious  approximate  approach 


to  determining  m^  and  m2  has  been  developed  and  will  be  dis¬ 
cussed  next.  The  approach  is  useful  in  a  transition  region  or 


any  region  where  superposition  of  two  or  more  current  compon¬ 


ents  results  in  some  curvature  of  the  experimental  In (I)  vs.  V 


curve . 
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AN  APPROXIMATE  APPROACH 


If  the  two-process  model,  Eq.  (16)  ,  is  valid,  and  if  the 
In (I)  vs.  V  curve  is  nonlinear,  one  may  obtain  a  pair  of  equa¬ 
tions  (24)  for  each  two  adjacent  data  points.  If  the  curve  is 
nearly  linear,  the  pair  of  equations  thus  obtained  may  not  be 
linearly  independent  or  may  not  have  a  solution,  so  that  the 
procedure  described  here  will  not  work  in  such  a  region.  If  a 
linearly  independent  pair  of  equations  is  obtained,  each  pair 
has  a  unique  solution  from  which  the  m's  can  be  determined  ap¬ 
proximately.  Since  each  such  pair  has  a  unique  solution 
whether  or  not  the  two-  or  one-process  model  is  valid,  values 
of  m  can  be  determined,  but  these  may  not  be  the  integer  values 
expected.  We  define 


G  = 


(32) 


RPM 


_1 _ 

m  ^2 


(the  reciprocal  product  of  m's) ,  and 


SRM 


m^  +  m2 
mlm2 


(33) 


(34) 


(the  sum  of  reciprocals  of  m's).  These  symbols  are  useful  com¬ 
puter  words.  Let  I  and  I  +  1  be  adjacent  data  points  and  look 
for  the  Jth  solution  of  the  pair  of  equations  obtained  from 
Eq.  (24).  To  be  compatible  with  computer  language,  let  Bf  -♦  BF 
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and  (gfj  ■*  BP2.  Also  rewrite  Eq.  (24)  as 


G  =  1  - 


(m^  +  m, 
mlm2 


-  Sf  -  (6f): 

m  ^2 


Then 


BF(I)SRM(J)  -  BF2 (I) RPM ( J)  =  G(I) 


BF (1+1) SRM ( J)  -  BF2 (1+1) RPM( J)  =  G(I+1) 


to  a  good  approximation.  The  solution  of  this  pair  of  equa¬ 
tions  is 


SRM(J)  = 


_  BF2  (14-1) G  (I)  - 


BF2  (1) G (1+1) 

bf(i+1)bf5 (I) 


RPM  ( J)  = 


_  BF (I+1)G (I)  -  BF (I) G (1+1) 


Letting  XM2  5  m2  and  XM1  =  m^. 


ymo/tn  _  SRM  ( J)  _  1 

XM2 2 RPM ( J)  2RPM(J) 


[SRMfJ);*  -  4RPM ( J)  (40) 


XM1(J)  =  HHtt  '  XM2(J)  * 


Once  the  m's  have  been  obtained,  there  are  nine  possible  combi 
nations  (which  will  not  be  enumerated)  of  results,  wherein  m. 
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1 


and  m2  are  constants  of  correct  value  (i.e.,  integers),  con¬ 
stants  of  incorrect  value,  or  variables.  If  one  of  the  m's  = 
the  current  component  may  be  injection  current.  If  so,  its 

correct  I  can  be  found  as  described  below.  If  the  other  m  is, 
o 

say,  2,  where  the  2  is  found  to  describe  space-charge  region 
recombination  current,  the  I  found  for  it  may  vary,  because 
I  for  such  current  is,  in  general,  injection-level  dependent. 
This  will  complicate  the  results  somewhat,  but  judicious  in¬ 
spection  may  help  in  deciphering  the  results .  For  any  case 

where  one  of  the  m's  is  not  constant,  the  value  of  I  obtained 

o 

may  be  treated  as  a  "subtotal"  current  of  the  form 


where  m  =  m(I) .  It  is  possible  to  treat  this  subtotal  current 
in  the  same  fashion  as  the  total  current:  Curve  fit  to  it,  as¬ 
sume  two  current  mechanisms  and  break  it  up  into  '  and  m2 ' , 
as  before,  continuing  until  the  currents  have  been  resolved 
satisfactorily.  The  foregoing  also  applies  to  any  case  where 
one  of  the  m's  is  constant  but  incorrect. 

In  any  case,  when  values  of  rt^  and  m2  have  been  obtained, 
we  write 

1  -  *01  exP<s7'  +  x02  exp'i7>  •  <43> 

For  two  adjacent  interpolated  data  points  J  and  J+l,  we  have 
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and  we  have  assumed  m^  and  m2  do  not  change  much  between  J  and 
J  +  1. 

To  demonstrate  two  of  the  approaches  which  have  been  used 
successfully^  and  to  test  their  validity,  data  from  a  hypotheti¬ 
cal  diode  were  analyzed  by  the  program.  This  diode  had  m^  =  1, 

“  8 

m2  =  2,  IQ1  =  10  units  and  IQ2  =  1  unit.  Its  I-V  character¬ 
istic  is  plotted  in  Fig.  1.  The  resultant  computer  printout 
is  given  in  Fig.  2. 

In  the  lowermost  and  uppermost  portions  of  the  curve, 
where  m2  and  m, ,  respectively,  dominate,  Eq.  (2)  is  an  excel¬ 
lent  approximation.  Hence,  XM  or  XMPRIM  describes  the  current 
mechanisms  very  well  and  one  could  deduce  the  components  easily 
for  this  diode.  PPRIM  is  equivalent  to  I  and  it  is  seen  to 
give  the  correct  values  at  either  end  of  the  range. 

For  the  transition  region  near  the  center  of  the  curve , 
the  value  of  m  as  calculated  from  Eq.  (2)  does  not  give  the 
correct  value,  as  expected.  In  this  region  (roughly  from  sub¬ 
scripts  I  =  50  to  I  =  70),  the  approach  using  Eqs.  (36)  through 
(41)  finds  its  application  and  the  values  of  m^  and  m2  are 
determined  with  good  accuracy  in  this  region.  Note  also  that 
Eq.  (27)  can  be  used  to  find  m^  and  m2  at  the  two  extremes, 
since  Sf  (called  BFP)  *  1  at  either  end  of  the  curve.  The  ex¬ 
perimental  results  from  real  diodes  usually  do  not  show  large 
regions  of  linearity  in  which  Eq.  (2)  is  useful,  so  that  the 
utility  of  the  latter  approach  in  regions  of  curvature  becomes 
very  apparent.  This  approach  does  not  work  well  in  the  linear 
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DIODE  VOLTAGE,  V,  volts 

Figure  1.  Current-voltage  Characteristic  of  a  Hypothetical 
Diode 


regions,  as  discussed  earlier.  The  values  obtained  in  the  non¬ 
linear  region  for  Iq^  and  Iq2  (called  CY01  and  CY02  or  CY01P 
and  CY02P  in  the  program)  are  not  as  satisfactory  as  those  ob¬ 
tained  for  m^  and  m2«  The  least  error  occurs  where  ir^  (or  m2) 
passes  through  its  correct  value.  In  this  case  the  error  may 
be  negligible  (less  than  1%) ,  but  it  may  be  off  by  a  factor  of 
two  or  three;  however,  it  is  well  within  an  order  of  magnitude 
in  the  worst  case,  and  this  would  often  be  adequate.  If  there 
is  a  linear  region  where  m^  or  m2  is  determined  correctly,  1^ 
and  Iq2  will  always  be  calculated  within  ±20%. 

At  present,  the  CAFA  program  is  in  a  primitive  stage  of 
development.  Further  work  should  provide  a  more  useful,  ac¬ 
curate  and  flexible  program  which  will  be  of  considerable 
value  in  experimental  research. 
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